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Abstract. We investigate bounds on the strength of the primordial magnetic field (PMF) 
from the cosmic microwave background (CMB) bispectra of the intensity (temperature) 
modes induced from the auto- and cross-correlated bispectra of the scalar and tensor com- 
ponents of the PMF anisotropic stress. At first, we construct a general formula for the CMB 
intensity and polarization bispectra from PMFs composed of any type of perturbation. Then 
we derive an approximate expression which traces the exact shape of the CMB bispectrum 
in order to reduce the computation time with respect to a large number of the multipole 
configurations, and also show that the non-Gaussian structure coming from PMFs is clas- 
sified as the local-type configuration. Computing the signal-to-noise ratio on the basis of 
the approximate formula with the information of the instrumental noises and resolutions, we 
find expected upper bounds on the magnetic field strength, when the magnetic spectrum is 
nearly scale invariant (hb = —2.9), smoothed on IMpc scale at 95% confidence level from the 
WMAP and PLANCK experiments as Bi Mpc < 4.0 - 6.7nG and 3.8 - 6.5nG, respectively, 
depending on the energy scale of the magnetic field production from 10 14 GeV to 10 3 GeV. 
Our new consequences imply slight overestimations by the previous rough discussions. 

Keywords: primordial magnetic fields, non-gaussianity, physics of the early universe, cos- 
mological perturbation theory 
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1 Introduction 

Thanks to the precise measurements of the cosmological and astrophysical events, there exist 
some inevitable evidences of the existence of the cosmological-scale magnetic fields, whose 
strength is about a few microgauss, in galaxies and clusters of galaxies at redshift z ~ 0.7—2.0 
[1-3]. In addition, recently, from the data of the Fermi and HESS experiments, it has been 
found that the magnetic field strength in the inter-galactic medium is larger than 0(1O~ 15 )G 
[4-6] . In order to realize the current magnetized Universe via the astrophysical magnification 
known for the dynamo mechanism (e.g., refs. [7-9]), a lot of theories predict the existence of 
the seed magnetic field in the early Universe (e.g., refs. [7, 10-12]). This seed field acts as 
a source of the cosmic microwave background (CMB) and generates characteristic patterns 
on its anisotropy. Therefore, the attempts to extract the properties of the seed field from 
the CMB two-point correlation function (power spectrum) have been well studied, and we 
had a broad understanding that the Universe favors the scale-invariant magnetic field whose 
magnitude (normalized by the present value) is less than 0(l)nG (see e.g., refs. [13-22]). 

On the other hand, even if the primordial magnetic field (PMF) obeys the Gaussianity, 
due to the quadratic dependence of the CMB fluctuation on them, the higher-order correlation 
of the CMB fluctuations also appear [23]. In particular, since the CMB three-point function 
(bispectrum) has been well utilized as a measure of the primordial non-Gaussianity (e.g., 
refs. [24-26]), there are a lot of literature which have been analyzing the impacts of the 
PMF on the CMB bispectrum and seeking its reasonable constraints [27-33]. However, in 
refs. [27-30], the authors approximately computed the integrals over the wave number vectors 
in the CMB bispectrum of the magnetic scalar modes. In our works [31-33], despite success 
of the exact calculation of the CMB auto-correlated bispectra between the vector and tensor 
modes, bounds on the PMF magnitude were estimated only by a rough comparison with the 
CMB bispectrum from curvature perturbations. In this sense, these estimations may have 
involved several kind of uncertainties. In addition, the previous studies have not focused on 
the contribution of the cross-correlated bispectra between the different perturbation modes 
yet, that is, the previous works seem to be incomplete. 
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The goal of this paper is to obtain bounds on the PMF strength by computing of the 
signal-to-noise ratio including all information of the scalar and tensor passive modes if the 
PMF spectrum has the nearly scale-invariant shape as ub = —2.9. Here, the passive modes 
denote the CMB signals induced by the growth of the metric perturbations due to the PMF 
anisotropic stress prior to neutrino decoupling, and dominate over the CMB fluctuation from 
PMFs for t < 1000 [15, 19]. To do this, we firstly give an exact formula for the CMB 
intensity and polarization bispectrum composed of any type of perturbation. However, as 
the numerical computation based on this formula requires a great deal of time, we present 
a faster form for the auto- and cross-correlated bispectrum between the scalar and tensor 
modes, which traces the shape of the exact one, by an appropriate approximation similar to 
the way of ref. [34]. Then, we also find that the bispectrum of the PMF anisotropic stresses 
represents the local-type non-Gaussianity. After that, using this approximate formula, we 
compute the signal-to-noise ratio in terms of the intensity-intensity-intensity bispectra with 
the effects of the instrumental noises and resolutions of the current observations as WMAP 
and PLANCK [25, 35], and obtain new bounds on the PMF strength. 

This paper is organized as follows. In the next section, we review the definition and 
statistical properties of the PMF. In section 3, we produce exact and approximate formulae 
for the CMB bispectrum, and analyze these behaviors. In section 4, we present the signal-to- 
noise ratio and the limits on the PMF strength. The final section is devoted to a summary 
and discussion of this paper. 

Throughout this paper, we obey the definition of the Fourier transformation as 



and the rule that the subscripts and superscripts of the Greek characters and alphabets run 
from to 3 and from 1 to 3, respectively. 



2 Statistics for the primordial magnetic fields 

In this section, according to ref. [32] , let us summarize the statistical properties of the stochas- 
tic primordial magnetic fields (PMFs), B(x, r), on the homogeneous background Universe 
characterized by the Friedmann-Lemaitre- Robertson- Walker metric as ds 2 = a(r) 2 (— dr 2 + 
Sb c dx b dx c ), where r is a conformal time and a(r) is a scale factor. The expansion of the 
Universe makes the amplitude of the magnetic fields decays as I /a 2 and hence we can draw 
off the time dependence as B(x, r) = B(x)/a 2 . Using this scaling relation, each component 
of the energy momentum tensor of the PMF is given by 




(1.1) 




8 b c -B\*)B c {*) =p 7 (r) A B (x*)# c 



S 2 (x) ee -p,(r)A B (xn 



+ n» Bc (x») 



(2.1) 
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where we have introduced the photon energy density p 7 to confine the time dependence of 
a -4 . After the Fourier transformation, the spatial parts are described as 

T 6 c (k,r) = p 7 (r) [5 b c A B (k) + U b Bc (k) 



As(k)= 8^y w™ k -^ ( 2 - 2 ) 

n fe Bc (k) = / -^^ (k / )Bc(k _ k / } ^ 



4vrp 7j0 7 (27r) 

where /9 7j o denotes the present energy density of photons. In the following discussion, for 
simplicity of calculation, we ignore the trivial time-dependence. Hence, the index can be 
lowered by 5b c and the summation is implied for repeated indices. 

Assuming that B(x) obeys the Gaussianity, its power spectrum is defined by 1 

(B a (k)B b (p)) = (2^) 3 ^lp a6 (k)5(k + p) . (2.3) 
Here, we use the projection tensor, which comes from the divergence free nature of PMFs, as 

P ab (k) = 4 CT) (k)4" CT) (k) = <U - kk , (2.4) 

£T=±1 

with k and being a unit vector and a normalized divergenceless polarization vector which 
satisfies the orthogonal condition; k a eit 1 ^ = 0. The details of the relations and conventions 
of the polarization vector are described in appendix A. Although the form of the power 
spectrum P B (k) strongly depends on the production mechanism of the PMF, we assume a 
simple power law shape given by 

P B (k) = A B k nB , (2.5) 

where Ab and n# denote the amplitude and spectral index of the power spectrum of PMFs, 
respectively. In order to parametrize the strength of the PMF, we smooth the magnetic fields 
with a conventional Gaussian filter on a comoving scale r: 

B 2 r = r^e-^P B (k). (2.6) 



Then, A B is calculated as 



(2vr) ns+5 B 2 r 
Y{^±l)k r r 



AB = -ns+S^ns+S > ( 2 ' 7 ) 



where T(x) is the Gamma function and k r = 2ir/r. 

As we will mention in the next section, the anisotropic stress of the PMF, H Ba b, acts as 
a source of the CMB fluctuation and hence we shall focus on its statistical properties. Using 
a projection tensor on the scalar mode 0^ b = — k a kb + ^afo/3, ei^, and a transverse-traceless 



1 Hore we neglect the helical term. 
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polarization tensor e ab = \f2e a ei (for more detail see appendix A), the fluctuation of 
the PMF anisotropic stress is decomposed into the scalar, vector and tensor parts: 



where in the second equality of the second equation, we use a fact that Usab is symmetric 
for the permutation of a and b. According to equation (2.1), the PMF anisotropic stress 
at an arbitrary point, H Bb (x.), depends quadratically on the Gaussian PMF at that point, 
B(x). This situation is identical to the local-type non-Gaussianity of curvature perturbations: 



C(x) = Cg(x) + (3/5)/^ al [C&(x) - <Cg( x )>] with CCg and /^ al being the curvature per- 



turbation, its Gaussian part and the nonlinear parameter of the local-type non-Gaussianity, 
respectively [36]. Hence, it is expected that the statistical properties of the PMF is similar 
to those of the local-type non-Gaussianity. This will be discussed in the next section. 

Using equation (2.3) and the Wick's theorem, the bispectrum of the PMF anisotropic 
stresses is calculated as 



where kr> is the Alfven-wave damping length scale [37, 38] as k~p ~ 0(O.l)Mpc and the 
curly bracket denotes the symmetric seven terms under the permutations of indices: a •<-)■ b, 
c <->■ d, or e o /. Here, for convenience of the calculation, we express as a more symmetric 
form than that of ref. [23]. To avoid the divergence of (UBabO<-i)^BcdO<-2)^BefO<-3)) in the 
IR limit, the spectral index is limited as tib > —3. Note that equation (2.9) is the six- 
point correlation of the Gaussian variables, that is, the higher-order correlation than the 
bispectrum of curvature perturbations. Hence, for formulating the CMB bispectrum induced 
by PMFs, appropriate treatments of the complicated computation corresponding to the one- 
loop diagram are required. 

3 CMB scalar, vector and tensor bispectrum 

In this section, we present the exact and approximate forms for the CMB bispectrum of the 
intensity and polarization modes induced by the scalar, vector and tensor perturbations of 
the non-Gaussian PMF anisotropic stresses, and discuss the consistency between the exact 
and approximate bispectra via numerical computations. 

Conventionally, the intensity, and E- and i?-mode polarizations (X = I, E, B) generated 
from the CMB scalar, vector and tensor anisotropies {Z = S, V, T) are expanded by the spin-0 
spherical harmonics as 




^(k^n^k) , 

\ [M Tl) ( k ) + ^4 Tl) (k)] n Ba6 (k) = ^e[ Tl) (k)n Ba& (k) , (2.8) 




,n=l 



n 



x<5(ki - k' x + k' 3 )5(k 2 - k' 2 + ki)<5(ks - 4 + k' 2 ) 

x^[P ad (k\)P be (K)Pcf(K) + b or c O d or e O /}], 

(2.9) 




(3.1) 
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where n is a line-of-sight direction. Then, the each-mode coefficient, a x I, is formally given 
by 

4\ m = M-tf / jffiS E[ s g n ( A )] A+ ^ (A) ( k ) r i?( fc )^^/ m (k) , (3.2) 

where A expresses the helicities of the scalar, vector and tensor perturbations as A = for 
Z = S, A = ±1 for Z = V and A = ±2 for Z = T, x discriminates the parities of the 
intensity and two polarizations as x = for X = J, E and x = 1 for X = B, and f W and 
T x J(k) are the primordial perturbation and transfer function of each mode, respectively. 2 
If PMF exists, its anisotropic stress acts as the source of the CMB anisotropy. For example, 
the primordial gravitational waves and curvature perturbations are amplified by the PMF 
anisotropic stresses prior to neutrino decoupling although these growths stop due to the 
compensation by the neutrino anisotropic stresses after neutrino decoupling. These extra 
gravitational waves and curvature perturbations survive passively, and generate additional 
fluctuations on the tensor- and scalar-mode CMB fields, respectively, similar to the non- 
magnetic ones on large and intermediate scales [15, 19, 39]. 3 On the other hand, in the 
vector mode, as the gravitational potential decays away via the Einstein equation, additional 
CMB fluctuation arises from not the above mechanism but the enhancement of the vorticity 
of photons by the Lorentz force from the PMF at small scales [13-15]. This Lorentz force is 
also originated from the PMF anisotropic stress, hence, we can express each-mode primordial 
perturbations with the PMF anisotropic stress: 

eW(k)*-iVn(^) ng(k), 

^(k^n^k), (3.3) 
e( ±2 >(k)«6ivn(-|) ngf(k). 

Here, note that and £( ±2 ) correspond to the superhorizon solutions of the Einstein 
equations with respect to the curvature perturbation and gravitational wave, 4 respectively 
[15, 19, 39, 42], and depend on the production time of the PMF, tb, the epoch of neutrino 
decoupling, t v ~ lMeV -1 , and the ratio by the energy density of photons to all relativistic 
particles, Ry ~ 0.6 for r < r v . In this paper, as the upper and lower values of r^ 1 , we take 
the energy scale of the grand unification (10 14 GeV) to the electroweak symmetry breaking 

(10 3 GeV), corresponding to t u /tb = 10 17 and 10 6 , respectively. As T X S J and 7^?, we can 
use the standard cosmological transfer functions without the magnetic dependence [43-46] 
because posterior to neutrino decoupling, the PMF contributes little to the evolution of cos- 
mological perturbations [19]. In contrast, as Txe' ', we should use the form involving the 
impacts of the PMF shown in Refs. [13, 14, 32]. 

2 Here, we set 0° = 1. 

3 In refs. [40, 41], by taking into account the effects in both the inflationary and the radiation-dominated 
epochs consistently, the authors concluded that the solution of the curvature perturbation spoils. See however 
ref. [42]. 

4 £(± 2 ) dgjjotgg eacn helicity state of the primordial gravitational wave, which is expressed as 

£ (±2) (k) = i e ^ 2) (k)Mk) , (3-4) 
with hij = Sgij/a 2 being the spatial components of the metric perturbation. 



- 5 - 



For convenience of the computation for the CMB bispectrum, we shall perform the 
multipole expansion of each-mode primordial perturbation as 



£W(k) = £^(A0-A^ m (k) . 

Im 



Then, equation (3.2) changes to 



i {Z) 

l Xlm 



(3.5) 



(3.6) 



and their bispectrum is easily formulated as 



if 

\n=l 



(Zu) 



n m-o* 



n=l 



oo i,2 



k n dk n ( Zn ) 



(2tt) c 



(3.7) 



\n=l 



In the following subsections, by expanding ^]^ =1 Ci^m„ (k n )^ and performing summations 
over Ai,A2, and A3, we firstly find an exact expression of the CMB bispectra from PMFs. 
However, it takes an awful lot of time to compute the CMB bispectra based on the exact 
form. Thus, using an approximate expression of the bispectrum of the PMF anisotropic 
stresses, we develop a faster form that can reconstruct the exact shape in section 3.2. 

3.1 Exact expression 

Here, we want to expand the CMB bispectrum (3.7) with the mathematically rigid method 
[32]. To obtain an analytic formula for the CMB bispectrum, we firstly should reduce the 
bispectrum of ^^{k), which is explicitly given by 



IT 

\n=l 



(An) ( 



n 

.71=1 



d 2 k n _ An y; nmn (k n ) 



n? {An) ( k n) 



(3.« 



\n=l 



Note that yHn=i ^ An ^(k n )y is composed of the bispectrum of the PMF anisotropic stresses 
(2.9) and the relations (2.8) and (3.3). To perform three angular integrals in this equation, 
we expand the delta function and projection tensors in /]^ =1 ^ A?l ^(k n )^) by using the spin- 
weighted spherical harmonics as 



\n=l 



y 2 dy 



xT 



000 



n e (-d w: 

n=l L n M„ 
L\ L2 L3 

Mi M 2 M 3 



Jl n M n ( k r 



(3.9) 



and some equations described in appendix A, respectively [32, 47]. Here, j\ is the Bessel 
function and the / symbol is defined by 



(2Zi + l)(2Z 2 + l)(2Z 3 + l) ( h h h 

4-7T I S i S 2 S3 



(3.10) 
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Furthermore, we replace the angular integrals of the spin spherical harmonics with the 
Wigner-3j symbols, and sum up these Wigner-3j symbols over the azimuthal quantum num- 
bers in the same manner as in ref. [32]. After these complicated treatments, we can obtain a 
formula for the initial bispectrum composed of any combination of three kind of the primor- 
dial perturbations as 



\n=l 



tl 12 C3 

m\ m2 rris 



(-47rp 7)0 



,-3 



n 

n=l 



k n dk n PB(k n ) 



ii i 2 «3 
L' L" L 



E E 

LL'L" S,S',S"=±1 

f S"SX 1(k > y k ,\fSS'X 2(k > ,/ h )fS'S"\ 3(k , 
JL"Lt x V 3' 1' K 1)J LL'i 2 \ K li K 2i^)J L'L"£ 3 l ft 2> ^3' ft 3 j 



(3.11) 



where 
Jl"U \rz,r2,r\) 



E 

L\L 2 l 3 

x(-r/ 



y 2 dyjL i (r 3 y)j L2 {r 2 y)jL 1 (ny) 



-L2+L3I 



■1) 



l i+ l 2+l 3 T os"-s" T os-s T ox-x 

1 LiL 2 L 3 1 L 3 lL" 1 L 2 IL 1 L X 12 




'_ 2 (Svr) 3 / 2 ^ In (r,/r B ) (A = 0) 



x < 



|(S) 3/2 A (A = ±1) 

^ -4(8^) 3 / 2 ,R 7 In {t u /t b ) (A = ±2) 



(3.12) 



Here, we stress that this expression involves all other modes than the vector- vector-vector and 
tensor-tensor-tensor ones which have been considered in our previous works [31-33]. Note 

( i\ h h \ 

that the confinement of mi , mo and m-? to the Winger-37 symbol ( ensures the 

rotational invariance of the CMB bispectrum, hence, in the following discussion, we focus on 
the angle- averaged form B^x^x^ \-ii 2 i 3 defined by 



IF 

\n=l 



IZn) 



h h h 

mi m2 m 3 



D 



(Z X Z 2 Z 3 ) 

x 1 X2X3,e 1 e 2 e 3 • 



(3.13) 



Substituting equation (3.11) into equation (3.7), we can obtain an explicit form of the angle- 
averaged bispectrum as 



tj{Z\Z 2 Z 3 ) 

n x 1 x 2 x 3 ,e 1 i 



n<- 

n=l 



2vr 2 



< 2 <P B «) 



LL'L" S,S',S"=±1 k J 
X /l,''lJi (^3) k'li ^i)/ll'^ 2 2 (^1^2 5 ^2) f yyi^ {k'2-.k^ k$) . 



(3.14) 



From this equation, we can find that the CMB bispectrum from PMFs is not obtained 
unless performing three summations with respect to the additional multipoles L,L', and L" , 
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tt it 




Figure 1. Diagrams with respect to the multipoles. The left panel corresponds to the multipolc 
configuration of equation (3.14). Due to the Wigner-6j symbol originated with the sextuplicate 
dependence on the Gaussian PMFs, l\,l<x and £3 arc linked with the dummy ones L, L', and L", and 
the one-loop structure is realized. The right panel represents the tree-structure diagram, which arises 
from the CMB bispectrum induced by the four-point function of the Gaussian seed fields. 



l\ 12 <-3 

which is correlated with ^1,^2 and l 3 via the Wigner-6j symbol < , „ >. These extra 

summations come from the sextuplicate dependence on the Gaussian PMFs. This situation, 
which corresponds to the one- loop diagram as shown in the left panel of figure 1, leads to a 
lengthy computation. In the next subsection, using an appropriate approximation, we reduce 
this loop calculation to the tree- level one (corresponding to the right diagram of figure 1). 

For the numerical calculation, we need to decompose the / functions in equation (3.14) 
and perform summations over S, S' , S" , A1 , A2 and A3. Using the analytic results as 

Ef 1 \L 2 +L'+L''+L 3 +L'+L" T 0S-S t 0S-S tOS'-S' tOS'-S' t 0S"-S" t 0S"-S" 
V L ) L's^L 1 L 2 1L L L'l\L> 1 L' 2 1L' 1 L 3 1L" 1 L%1L" 



S,S',S"=±1 



tOX-X rOl-1 rOl-1 rOl-1 rOl-1 rOl-1 

1 Tt 1 t 1 r.oi t. 1 r/'i T ' 1 T.i 1 r.' J 7",oi t." 1 t 



1 L' 3 XL 1 L 2 XL ± L'l XL' 1 L' 2 XL' ± L 3 XL" ± L% XL" 
x f 8 (L> 3 + L 2 ,L'> + L' 2 , L 3 + L'i = even) 
1 (otherwise) 



(3.15) 



Ai A2A3 



_ r 0|A 1 |-|A 1 | 7 -0|A 2 |-|A 2 | r 0|A 3 |-|A 3 | 
_ 1 L 1 h 2 I L[ i 2 2 1 L'{ £ 3 2 

x f 2 3 ~ N s (Li + t x + Xl ,L[ + £ 2 + x 2 , L'{ + £ 3 + x 3 = even) 
I (otherwise) 

with Ns being the number of the scalar modes constituting the CMB bispectrum, 5 equation 
(3.14) is rewritten as 

BxixISaMz = c z x Cz 2 Cz 3 (-4tt / 9 7i0 )" 3 { l' L" L } 

T.TJ T .'I ^ * 



£ 

L\L 2 L 3 
L i L 2 L 3 

L 'l L 2 L 3 



LL'L 

L n +L'+L"+2t 



I 1 \E„-i " 9 " rOOO rOOO rOOO 

1 L 1 L 2 L 3 1 L' 1 L' 2 L' 3 1 L'{L%L% 



'Caution about a fact that |A| is determined by Z, namely, |A| = 0, 1,2 for Z — S,V,T, respectively. 
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K n aK n n -(Z n ) /, \ 
2yr 2 'X n ,£ n ^ Kn ) 



A dAjLi {k\A) / B A dBj L ,(k 2 B) 



C 2 dCj L ,,(k 3 C) 



k^dkiPBik'^jLMA^L'^B) I k> 2 2 dk' 2 P B (k> 2 )j L ,(k' 2 B)j L ,,(k' 2 C) 

k' 3 2 dk' 3 p B (k' 3 )j L ,Ak' 3 c)j L MA) 



x2 



oi-i 7-01-1 7-01-1 T-oi-i roi-i T-oi-i n n n 

L' 3 lL 1 L 2 lL 1 L'.llL' 1 L' 2 lL< 1 L i lL" 1 L>>lL" ^L' 3 ,L 2 ,L^L' 3 \L' 2 ,L'^L 3 ,L^L" 

(e) 



3-7V S7 -0|A 1 |-|A 1 |.0|A 2 |-|A 2 | r 0|A 3 |-|A 3 | r) ( e ) ^(e) ^ 



L[ £i 2 

Here, we define the filter functions as 



V,X\ L[,£2,X2 L'J,£3,X3 



Ql' 3 ,l 2 ,l = (8l' 3 ,l+i + ^,|l-i|)(<5l2,l+i + <5l 2 ,|l-i|) + 5l' 3 ,l5l 2 ,l 
which come from the selection rules of the Wigner symbols [47], and 



(3.16) 



(3.17) 



Cz 



--|(8tt) 3 / 2 P 7 In (t„/t b ) (Z = S) 
|(8tt) 3 / 2 



(Z = V) 



(3.18) 



-4(8vr) 3 / 2 P 7 In (t u /t b ) (Z = T) 
This expression is used in the calculations of the exact bispectra presented in section 3.3. 
3.2 Pole approximation 

As we have mentioned in the previous subsection, it takes a lot of time to compute the CMB 
bispectrum based on the exact formula corresponding to the one-loop calculation. Therefore, 
in order to constrain on the PMFs through the estimation of the signal-to-noise ratio in 
a finite time, we find the reasonable expression of the CMB bispectrum by reducing the 
exact formula to the tree-level one. A similar approach can be seen in the calculation of the 
non-Gaussianity generated from scalar-type Gaussian-squared fields such as curvatons [34]. 

Let us consider the structure of the bispectrum of the PMF anisotropic stresses (2.9). 
This consists of three integrals over the dummy wave number vectors k'-^k^, and kg. If the 
magnetic spectrum is enough red as n B ~ —3, these integrals are almost determined by the 
behavior of the integrand around at "three poles", namely, ki,k 2 ,k 3 ~ 0. Picking up the 
effects around at these poles, we can express the bispectrum of the PMF anisotropic stresses 
with an asymptotic form as 



(n Safe (ki)n Bcd (k 2 )n Be/ (k 3 )) ~ (-47rp 7i0 ) 



-3 aA B k n B +3^ § 
n B + 3 * 3 



\n=l / 



l r 



p B (fe 1 )p B (A 2 )J ad p 6e (k 1 )p c/ (k 2 ) 
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+P B (k 2 )P B (k 3 )P ad (k 2 )P be (k 3 )6 cf 
+P B (k 1 )P B (k 3 )P ad (k 1 )6 be P cf (k 3 ) 
+{a oiorc«(ioref)>/}] , 



where we have evaluated as 



fe» 



d 3 k'P B (k')P ab (k>) ~ a / k' 2 dk'P B (k') / d 2 k'P a6 (k') 



(3.19) 



(3.20) 



Here, we introduce k* = lOMpc" 1 and a as a cutoff scale of the integral and an unknown 
numerical factor which should be determined by the comparison with the exact bispectra 
described in the previous section, respectively 6 . 

Once we admit this validity, we formulate the four CMB bispectra generated from the 
scalar and tensor components of the PMF anisotropic stresses, namely the SSS, SST, STT, 
and TTT spectra 7 . The angular bispectrum of the primordial perturbations is approximately 
given by 
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(Ai = A 2 = 0, A 3 = ±2) 



T°i?( k >cd A2) (k^7 A3 W (Ai = 0, A 2 , A 3 = ±2) 
-27ei; Al) (k 1 )e(- A2) (k 2 )e(7 A3) (k 3 ) (A lf A 2 , A 3 = ±2) 



(3.21) 



Using some relations in appendix A, we reduce the contraction of the subscripts in the SSS 
spectrum to 



°ab (ki)o£? fo)0f) (k 3 ) \P B (ki)P B (k 2 )5 ad P be (ki)P c/ (k 2 ) 
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£ £ £ I t01-1 t01-1 r01-l 

1111 in L ' n L " n 



'A/'( k 2)^ / L»j\/»(k 3 ) 



MM'M" 



L L' L" 
M M' M" 



6 Equation (3.20) is valid for ub 7^ —3. 

7 Of course, we also can apply the pole approximation to the CMB bispectra composed of the vector modes. 
However, we do not focus on these bispectra here since the contribution of the vector mode is subdominant 
for the scale where the current satellites are sensitive. 
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x [PB(h)PB(k2)8 L n !2 + SL, 2 P B (k 2 )P B (k 3 )+PB(k 1 )S L , t2 P B (k 3 )] , (3.22) 
that for the SST spectrum to 

( " b \ki)0^(k 2 )e { - X3 \k 3 ) [p B (h)P B (k 2 )5 ad P be (ki)P cf (k 2 ) 
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M M' M" 



X 

MM'M" 

x [P B (fci)PB(fc2) + 3{5 l , 2 Pb(^) + PB(h)S L/>2 }P B (k 3 )] 5 L „, 2 , (3.23) 



that for the STT spectrum to 

0^(k\)e[-/ 2 \k 2 )e { - X:i \k 3 ) [Pu(A:i)P B (fc 2 )^P f)e (k 1 )P c/ (k 2 



/ 
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L L' L" 
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x 

MM'M" 

x [P B (A;i){P s (fe 2 ) + P B (k 3 )} + 3<5 Li2 P B (fc 2 )P iJ (A;3)]<5 L , i2 <5 L » i2 , (3.24) 



and that for the TTT spectrum to 



'(k 1 )e^ A2J (k 2 )e^ A ^(k 3 ) [P B (k 1 )P B (k 2 )5 ad P be (^i)Pcf(k2) 

+Ps(A: 2 )Pi 3 (A:3)P ad (k 2 )P 6e (k3)5 c/ + PB(A; 1 )P B (A ; 3)P ad (k 1 ) ( 5 be P c/ (k3; 

\ / t T i Tii—n i y ' 



L,L',L"=0,2 

(L V L" 
M M' M" 

1V1,1V1 ,1V1 

x[P B (ki)P B (k 2 ) + P B (k 2 )P B (k 3 ) + P B {k 1 )P B (k 3 )}8 L , 2 8 va 8 L „2 . (3.25) 

The delta function is also expanded with the spin spherical harmonics as equation (3.9). 
Then, the angular integrals of these spin spherical harmonics are performed as 

f fik, , Y* Y* ^ Y* - r A i°~ A i ( ^ Ll L \ 
J a 1— Ai "iimi ^LiA/i Ai * LM ~ h x L x L \ mi M X Mj ' 

d2 ^-M Y t 2 m 2 Y L 2 M 2 \2 Y L' M' = ^£ 2 L 2 L' 2 ]tf 2 M' ) ' ^' 26 ' ) 

W 2 lA , V* V* , V* — 7~ A 30-A 3 / ^3 ^3 L 
a K-Z-X-s 1 Ism?, 1 LzMzM 1 L" M" ~ 1 " 



i 3 L 3 L" \ m3 M3 M u 
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and resulting Wigner-3j symbols are summed up as 
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Thus, the resultant initial bispectrum (3.21) are rewritten as 




(3.27) 
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(3.28) 



Substituting the above expression into equation (3.7) and using the result of the summation 
over Ai, A 2 and A3 (3.15), the resultant approximate CMB bispectrum composed of the scalar 
and tensor modes is derived: 
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y 2 dy 
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x < 



3[P fl (A: 1 )Pfl(fe2)ffi//,2 
+<5l, 2 Pb(A:2)Pb(A:3) 
+PB(fci)<5zy, 2 PB(fc 3 )] 

6V3[P B (fc 1 )P i? (A;2) 

+3{*L j2 Pfl(fe) + P B (A; 1 )5 L / i 2} 
xP B (fc 3 )]5 L // i2 

108[P B (A; 1 ){P B (A ;2 ) + P B (A ; 3)} 
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(3.29) 



Comparing the above expression with the exact form of the CMB bispectrum [(3.14) or 
(3.16)], we can see that the numbers of the time integrals and dummy multipoles decrease. 
In addition, note that the dependence on ki,k 2 and k 3 is identical to that in the CMB 
bispectrum arising from the local-type non-Gaussianity of curvature perturbations which is 
given by 



B 



local 



rOOO 
1 tilth 



y 2 dy 



n 

.n=l 



k^dknT^ j (k n )ji n (k n y) 



xhhT l [Pdki)Pd 



+ 2 perms.] 



(3.30) 



with /nl 3 " 1 an d Pf being the so-called nonlinearity parameter of the local-type non-Gaussianity 
and the power spectrum of curvature perturbations, respectively. These imply that by the 
pole approximation, the one- loop calculation (the left panel of figure 1) changes to the tree- 
level one (the right one of figure 1). Thus, we can reduce the time cost involved with 
calculating the bispectrum. 

3.3 Numerical issues 

Here, we present the numerical results of the exact and approximate bispectra for n B = —2.9 
based on equations (3.16) and (3.29). To compute the actual bispectra, we modified the 
Boltzmann Code for Anisotropics in the Microwave Background [15, 48] and use the Common 
Mathematical Library SLATEC [49]. 

In figure 2, we plot the reduced bispectra of the intensity modes, which is given by 



, {Z 1 Z 2 Z 3 ) = (T o \-l MZiZ 2 Z 3 ) 

J ni,he 2 e3 — ^hhh) D in,t 1 i 2 l ? , ' 



(3.31) 



on the basis of the exact (circle points) and approximate (lines) formulae of the magnetic 
case for (Zx, Z 2 , Z 3 ) = (S, S, S), (S, S, T), (S, T, T), (T, T, T), (V, V, V), and the non-magnetic 
formula (3.30). From this figure, we can see that the approximate SSS, SST, STT and TTT 
spectra (3.29) are in good agreement with the exact ones (3.16). Quantitatively, using the 
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Figure 2. Absolute values of the CMB bispectra of the intensity fluctuations generated from the 
magnetic SSS, SST, STT, TTT and VVV modes, and the local-type non-Gaussianity of curvature 
perturbations (3.30) with /^l = ^ f° r £1 = £2 = £3 = £■ The circle points and lines of the 
magnetic SSS, SST, STT,TTT and VVV modes correspond to the exact bispectra (3.16) and the 
approximate ones (3.29), respectively. The PMF parameters are fixed to -BiMpc = 4.7nG, tib = —2.9 
and t v /tb = 10 17 , and other parameters are identical to the mean values derived from the WMAP 
7-yr data [25] 



correlation function for each mode b • b' = YLt^ii^m^^n^Ul^ •> ^ ne correlation coefficients 
between the exact (b ex ) and approximate (6 app ) bispectra are evaluated as 
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(3.32) 



These become almost one, therefore in the following discussion, we believe that equation (3.29) 
reconstructs the shapes of the exact bispectra. 

The remaining factors, a's, which are introduced in equation (3.20), are determined by 
an asymptotic relation as 
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(3.33) 



These values are used in figure 2 and in the estimation of the signal-to-noise ratio described 
in the next section. 

From here, let us discuss the behaviors of these bispectra. On large scales (£ < 100), the 
tensor mode is the most significant signal in the intensity anisotropy due to the boost of the 
Integrated Sachs- Wolfe effect [33, 50], and highly damped at £ > 100. On the other hand, the 
scalar mode stands out at I > 100 owing to the acoustic oscillation of the baryon-photon fluid. 
Although the damping by the viscosity of photons is effective for I > 1000, the scalar mode 
keeps on dominating the tensor one [19]. The vector mode is monotonically increasing for 
£ < 2000 due to the vorticity of photons enhanced by the Lorentz force [13, 14]. Regardless of 
it, this acts as a subdominant signal for £ < 1000 [19]. These features can be seen in the CMB 
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bispectra. That is, from figure 2, it is observed that the more tensor modes the bispectrum 
comes from, the larger scale it shows up on. In this figure, the contribution of the vector 
mode seems to be negligible. In addition, we also notice that the magnetic SSS bispectrum 
and the bispectrum from the local-type non-Gaussianity of curvature perturbations have the 
similar shapes with each other despite the different dependence on multipoles. This backs 
up the fact that the non-Gaussianity of the PMF anisotropic stresses is classified into the 
local-type configuration. 

4 Signal-to-noise ratio 

In this section, using an approximate formula for the CMB bispectrum generated from the 
scalar and tensor perturbations of the PMF anisotropic stresses (3.29), we compute the signal- 
to-noise ratio and estimate the optimal observational bound on the PMF strength when the 
PMF spectrum is nearly scale invariant as ub = —2.9. 

According to refs. [24, 36, 51], the signal-to-noise ratio with respect to the CMB bis- 
pectra is formulated as 



where A^ 2 ^ 3 = 1 (for 4 / ^ / ^3), 6 (for l\ = £2 = £3), 2 (otherwise) denotes a numerical 
factor caused by the cosmic variance, Cf d is the fiducial CMB power spectrum which is 
consistent with the current observational data [25], and iVg is the noise spectrum given by 



with #fwhm and AT being the full width at half-maximum of the Gaussian beam and the 
instrumental noise per pixel in terms of temperature, respectively. 

Firstly, as an ideal limit, we describe the signal-to- noise ratios derived from the SSS 
(green line), SST + STS + TSS (magenta one), STT + TST + TTS (blue one) and TTT (red 
one) spectra and their total spectrum (orange one) when N( = in figure 3. From these, we 
can analyze the behavior of the total signal-to-noise ratio: at first, owing to the boost of the 
TTT spectrum, the total signal-to-noise ratio rapidly increases for l max < 100. Below that 
angular scale, the contribution of the TTT spectrum is completely saturated for ^ max ^ 100, 
and the cross-correlated bispectra between the scalar and tensor modes lead to the gradual 
growth. If £ m ax reaches 700, the contribution of the SSS spectrum becomes effective and the 
total signal-to-noise ratio starts to increase again. One may wonder about a fact that S/N 
from the total CMB spectrum is smaller than that from the TTT spectrum at l ma , x ~ 100. 
This is due to the compensation of the total CMB bispectrum induced by the other-mode 



Figure 4 shows the signal-to-noise ratios if we take into account the finite instrumental 
noises and angular resolutions in the WMAP and PLANCK experiments, respectively as 
(#fwhm [arcmin] , AT[fiK]) = (13, 18.3), (5.0, 13.1) [53-55]. We can see that the signal-to- 
noise ratios with the WMAP and PLANCK noises start to peak out from I ~ 500 and ~ 1000, 
respectively. From this figure, we shall find the observational limits on the PMF magnitude for 




(4.1) 



[51, 52] 




(4.2) 



spectra. 
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Figure 3. Noise-free signal-to-noise ratios. The "total" line denotes S/N obtained from the total 
spectrum of the SSS, SST, STS, TSS, STT, TST, TTS and TTT spectra, and the others correspond 
to S/N's coming from each mode. Here, we fix the PMF parameters as Ub = —2.9. The other 
parameters are identical to the mean values obtained from the WMAP-7yr data [25]. 
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Figure 4. Signal-to-noise ratios with zero, PLANCK- and WMAP-levcl noises obtained from the 
auto- and cross-correlated bispectra between the scalar and tensor modes. The instrumental noises 
and resolutions of the WMAP and PLANCK experiments are given by (#fwhm [arcmin] , AT[/iK]) = 
(13, 18.3), (5.0, 13.1), respectively [53-55]. Here, we fix the PMF parameters to the values presented 
in figure 3 and the other parameters are identical to the mean values obtained from the WMAP-7yr 
data [25]. 



n B = -2.9. At Anax = 1000, the signal-to-noise ratios with the WMAP and PLANCK noises, 
which are normalized by -BiMpc — InG and t u /tb = 10 , are 5.13 x 10~ 4 and 6.16 x 10~ 4 , 
respectively. Then, equating these values with S/N oc (-BiMpc) 6 '[^(tu/tb)] 3 , we can see that 
S/N < 2 and t v /t b = 10 17 - 10 6 lead to the upper bounds on the PMF strength at 95% 
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confidence level: B 1Mpc < 4.0 - 6.7nG (WMAP limit), 3.8 - 6.5nG (PLANCK limit). One 
may be disappointed that these bounds are somewhat relaxed compared with the previous 
ones estimated from only the SSS or TTT spectrum [29, 33]. However, these were based 
on a rough estimation and may involve some uncertainties, hence, our new results are more 
reliable. 

5 Summary and discussion 

In this paper, we newly studied how the PMF strength is constrained from the auto- and cross- 
correlated bispectra between the CMB intensity fluctuations of the scalar and tensor passive 
modes in the case where the PMF spectrum closes to the scale invariance as = —2.9. To 
do this, at first, we derived the most general formula for the CMB bispectrum generated from 
the bispectrum of the PMF anisotropic stresses including all considerable correlations with no 
approximation in accordance with the computation approach of ref. [32] . However, due to the 
dependence of the bispectrum of the PMF anisotropic stresses on the six-point correlation of 
the Gaussian PMFs, this formula involves a huge computational effort corresponding to the 
one-loop diagram. To avoid this difficulty, by evaluating the behaviors of the integrands with 
respect to the dummy wave number vectors around the poles, we approximately expressed the 
bispectrum of the PMF anisotropic stresses with the four-point correlation of the Gaussian 
variables and reduced the formula to the tree- level bispectrum. In this attempt, we also 
confirmed that the non-Gaussianity originated from the PMF anisotropic stresses has the 
local-type configuration. Through the comparison with the exact spectra, we checked that 
this approximate formula appropriately produces the exact shape. 

Using the approximate formula, we estimated the signal-to-noise ratio. For the noise- 
free case, we analyzed the behavior of this and checked the consistency of our calculation. 
Considering the instrumental noise and angular resolutions of the WMAP experiment, we 
obtained a bound on the PMF strength at 95% confidence level: -Bim P c < 4.0 — 6.7nG, 
where the upper and lower values correspond to the cases of t u /t b = 10 17 ,10 6 . In the 
PLANCK experiment, this value becomes slightly tigher as -BiMpc < 3.8 — 6.5nG. One can 
notice that these limits are weaker than the recent constraint and PLANCK forecast by the 
CMB power spectra reported in ref. [22]. The analysis in ref. [22] also includes information 
about the polarizations and contributions of the vector mode, which we do not consider 
here. Constraining from the polarization bispectrum induced by PMFs will be an interesting 
subject and one of our future works. 

In the previous attempts about probing the PMF with the CMB bispectrum [27-29], 
there exist many uncertainties in the evaluation of the CMB bispectra. In our recent works 
[31-33], regardless of the precise treatments of the CMB bispectra, the bounds have been 
roughly estimated. In this sense, we note that this is the first paper in which the properties 
of the PMF are extracted from the CMB bispectrum in the appropriate way. This study will 
be extended to the higher-order correlations such as the trispectrum [56]. If this realizes, we 
will acquire deeper comprehensions of the nature of the PMF. 
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A Projection tensors 

Here, we briefly summarize the useful properties of the projection tensors of the scalar, vector 
and tensor modes on the basis of ref. [47]. 

At first, we define an arbitrary unit vector, polarization vector and tensor, which are 
expanded by the spin spherical harmonics, as 



m 

6i ±1 )(k) = T ^c ±1 y lm (k) , (A.i) 



e it 2) (k) = V2 e (^(k) e r(k), 

with 

-m(6 m) x + 8 m -i) 

m,l + 

The scalar product of a™, which is required for the tensor contraction, is easily calculated as 



a T = \lY I i 0W + <W-i) I • (A.2) 



a a a a = "§"( _1 ) °m,-m' , = ^O^m' ■ (A.3J 

These polarization vector and tensor obey the divergenceless and transverse-traceless 
conditions and some properties in the Fourier space as 

^ei ±1 )(k) = 0, 

4 ±1) *(k) = ^ 1) (k) = e i ±1) (-k), 

4 A) (k)6i A ')(k) = 5 A ,_ A , (for A,A' = ±1) , 

eit 2) (k) = Uf ) (k) = 0, (A ' 4) 

4f ) *(k) = e if) ( k) = eif)(-k), 

e^(k)e^(k) = 26x,- x , (for A, A' = ±2) . 

Using these vectors and tensor, an arbitrary physical tensor, such as the metric or the PMF 
anisotropic stress, is decomposed into the two scalar (xiso,X^), two vector (x^ 1 ^) and two 
tensor (x^ 2 ^) components: 



Xab (k) = -- Xiso (k)5 a6 + x W(k)OW(k) 

+ £ X (A) (k) [k a e { b X \k) + k b ei X \k)] + £ X W(k)e^(k) , (A.5) 

A=±l A=±2 



- 18 - 



where ® h (k) = —k a kb + 5 ab /3. Considering equation (A. 4), we have the inverse formulae as 
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(A.6) 



From these, we find several relations of the projection tensors as 
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These representations are useful for the tensor contractions under the exact and pole approx- 
imation as shown in sections 3.1 and 3.2. 
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